Effective zero-thickness model for a conductive membrane driven by an electric field 
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The behavior of a conductive membrane in a static (DC) electric field is investigated theoretically. 
An effective zero-thickness model is constructed based on a Robin-type boundary condition for the 
electric potential at the membrane, originally developed for electrochemical systems. Within such 
a framework, corrections to the elastic moduli of the membrane are obtained, which arise from 
charge accumulation in the Debye layers due to capacitive effects and electric currents through the 
membrane and can lead to an undulation instability of the membrane. The fluid flow surrounding 
the membrane is also calculated, which clarifies issues regarding these flows sharing many similarities 
with flows produced by induced charge electro-osmosis (1CEO). Non-equilibrium steady states of 
the membrane and of the fluid can be effectively described by this method. It is both simpler, 
due to the zero thickness approximation which is widely used in the literature on fluid membranes, 
and more general than previous approaches. The predictions of this model are compared to recent 
experiments on supported membranes in an electric field. 

PACS numbers: 87.16.-b, 82.39.Wj, 05.70.Np 



I. INTRODUCTION 

Bilayer membranes formed from phospholipid 
molecules are an essential component of the membranes 
of cells. The mechanical properties of equilibrium 
membranes are characterized by two elastic moduli, the 
surface tension and the curvature modulus which 
typically depend on the electrostatic properties of the 
membranes [2|. Understanding how these properties 
are modified when the membrane is driven out of 
equilibrium is a problem of considerable importance to 
the physics of living cells. A membrane can be driven 
out of equilibrium in many ways, for instance by ion 
concentration gradients or electric fields, either applied 
externally or generated internally. 

The external application of electric fields on lipid films 
is used to produce artificial vesicles (by electroformation) , 
as well as to create holes in the membrane (by electro- 
poration) [3J. Both processes are important for biotech- 
nological applications, they are widely used experimen- 
tally although they are still rather poorly understood. 
The generation of ion concentration gradients by inter- 
nal means is controlled in biological cells by membrane- 
bound ion pumps and channels, which play key roles in 
many areas of biology [J] . 

The non-equilibrium fluctuations of membranes in- 
cluding ion channels and pumps were first analyzed in 
Refs. [1, H| by means of an hydrodynamic theory. Arti- 
ficially made active membranes inspired by these ideas 
were then studied experimentally 0, H, . Several the- 
oretical studies followed, mainly motivated by the ques- 
tion of how to model non-equilibrium effects produced by 
protein conformation changes [lfj, [llj, |l2j, ll3| . One lim- 
itation of existing active membrane models is that they 
do not describe electrostatic effects associated with ion 



transport in details. In previous papers by our group 
[3, [l5[ , we have addressed this limitation by studying a 
theoretical model for a membrane with a finite conduc- 
tivity transverse to the membrane plane (due for instance 
to ion channels or pumps) using electrokinetic equations 
[Trj . Il7l . [l8j . Our work complements Ref. (HI, where 
the correction to the elastic moduli of a membrane in 
a DC electric field were calculated using an approach 
purely based on electrostatics (no currents). It is also 
inspired by Ref. (2(| and Ref. [2l[ where similar prob- 
lems were considered using electrokinetic equations. In 
contrast to these studies, our approach focuses on the 
non-equilibrium case, where electrokinetic corrections to 
the elastic moduli arise due to currents through the mem- 
brane. 

In particular, a negative correction to the surface ten- 
sion arises due to capacitive effects, also called Lippman 
tension [22]. This negative tension leads to instabilities 
as can be understood from the high-salt limit (23|. A 
first experimental proof of the destabilizing effect of the 
electric field on a stack of lipid membranes was brought 
by X-ray scattering studies [Hj|. Recently, the lowering 
of the tension due to electrostatic or electrokinetic effects 
has been observed experimentally with supported mem- 
branes subjected to an AC electric field [25| and in active 
membranes @. 

The resulting flow fields around the undulating mem- 
brane are interpreted within the framework of 'induced 
charge electro-osmosis' (ICEO) [IE SI, [27j] . Similar flow 
patterns within vesicles subject to AC electric fields have 
been observed experimentally and analyzed theoretically 
in Ref. [28[ . The deformation of lipid vesicles in alternat- 
ing fields in various medium conditions has been modeled 
theoretically in Refs. [29|,|30(. All these studies show that 
lipid membranes in electric fields present a rich panel of 
possible behaviors 0, [3l| . 
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This paper extends previous work [lj, yj5( , by provid- 
ing an effective zero-thickness membrane model that con- 
tains both capacitive effects and ionic currents. In a first 
attempt [Tij , a zero-thickness membrane model has been 
proposed with the boundary condition (BC) of zero elec- 
tric field at the membrane. Although the shape of the 
potential was acceptable, the charge distribution had the 
wrong sign and the elastic moduli were orders of magni- 
tude too small. A model with finite membrane thickness 
and dielectric constant has thus been considered in Ref. 
fl5l ]. leading to correct signs of the charge distribution 
and correct orders of magnitude of the elastic moduli. 
However, this model needed approximations and finally 
numerical evaluations. In view of this, we present here 
an improved zero-thickness model, by using the more re- 
alistic BC of a dielectric interface sustaining Faradaic 
currents [32|. Although this Robin- type BC has been 
introduced in Ref. [15| , its consequences were not devel- 
oped. In particular this model leads to simple analytical 
expressions for the corrections to the elastic constants of 
the membrane. The model clearly captures both non- 
equilibrium effects due to ion currents and equilibrium 
effects, of capacitive nature. We also calculate the flow 
field around the membrane, which has in fact the oppo- 
site sign as compared to the one of Ref. [l5[ for the zero- 
thickness case, and is thus similar to standard ICEO flow 
fields. The presented effective zero-thickness model for a 
DC-field driven conductive membrane is simple enough 
to be the starting point of more refined further studies. 

The work is organized as follows: in section [TT] we de- 
scribe the equations for the charges in the electrolyte. 
A special emphasis is put on the boundary conditions 
which is the crucial point here. Then the base state so- 
lution corresponding to a flat membrane is calculated in 
section Mil In section IIVI we calculate the leading or- 
der contributions to the electric and ion density fields 
for a spatially modulated membrane height and analyze 
the corresponding hydrodynamic flows around the mem- 
brane. Using the boundary conditions for the stress ten- 
sor at the membrane (which includes Maxwell and hydro- 
dynamic stresses) , we calculate in section [V] the growth 
rate of membrane fluctuations. In section |VH our results 
are discussed and compared to previous calculations and 
to related experiments. 



distributions of channels / pumps in membranes we refer 
to Refs. [Ill, [Hj]. A point in the membrane is charac- 
terized by its Monge representation (valid in the limit 
of small undulations) by introducing a height function 
h{Tj_), where is a two-dimensional in-plane vector. 

In the electrolyte, the governing equation for the elec- 
tric potential <fi is Poisson's equation 
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Here e is the elementary charge and e is the dielectric 
constant of the electrolyte. For the sake of simplicity, 
we assume a symmetric 1 : 1 electrolyte, so that far away 



from the membrane n" 1 



We also assume that 



the total system is electrically neutral. 

The densities of the ion species are assumed to obey 
the Poisson-Nernst-Planck equations for a dilute solution 
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with ionic current densities 
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where fcgT is the thermal energy. We have assumed 
that both ion types have the same diffusion coefficient 
D, and neglected various corrections for concentrated so- 
lutions [3J]. We consider a steady state situation and 
use the Debye-Huckel approximation by linearizing the 
concentrations rv^ — n* + 5n^ , leading to 
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For symmetric binary electrolytes, it is useful to intro- 
duce half of the charge density [32|, [35( , 
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II. MODEL EQUATIONS 

We consider a steady (DC) current driven by a voltage 
V between two electrodes at a fixed distance L, applied 
to an initially flat membrane located initially at z = 0. 
The membrane is embedded in an electrolyte of mono- 
valent ions with densities n + and n~ . The membrane 
has channels for both ion species but is itself neutral (no 
fixed charges at the membrane). The channels/pumps 
are assumed to be homogeneously distributed in the 
membrane and enter only in the effective conductance 
G, as introduced below. For the effect of non-uniform 



as well as the average concentration of the two ionic 
species, 



5n^ 



Sn 



The latter quantity turns out not to be a relevant vari- 
able in the following since it is decoupled from the field at 
small applied voltages [35|. Moreover, since we have con- 
sidered a steady state and a symmetric situation, there 
is no net particle current. We also define 
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which represents half of the electric current density, and 
we arrive at the equations 
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Insertion of Eq. (fTO]) into Eq. (fTTjl leads to 

(d 2 z - k 2 ) p = , (12) 

where 

2e 2 n* 



ek B T 



(13) 



and k _1 = Ad is the Debye length that defines the char- 
acteristic length scale for charge relaxation in the elec- 
trolyte. 



A. Boundary conditions 



At the electrodes located at z 
impose the voltage leading to 



±-&, we externally 
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(14) 



This BC is oversimplified for real electrodes, since it 
neglects interfacial polarization across the double layers 
passing Faradaic currents 32, 36], which makes the volt- 
age imposed across the electrolyte, outside the double 
layers, different from the applied voltage. Since we focus 
on the membrane dynamics, however, electrode polariza- 
tion is inconsequential, and the voltage V in the model 
simply serves as a means to apply a steady DC current, 
which could be directly measured or imposed in exper- 
iments testing our theory. We assume in the following 
that the distance between the electrodes is much larger 
then the Debye length, L 3> Ad = k . In that case, the 
bulk electrolyte is quasi-neutral, n + = n~~ = n* , with 
negligible charge density (compared to the total salt con- 
centration) , 



P z 
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L 



= 0. 



(15) 



Since the conductivity of a quasi-neutral electrolyte is 
constant, the applied uniform current is equivalent to an 
applied electric field far from the membrane. 

As we will see, the BC at the membrane is crucial to 
recover the correct physical behavior. In the simple zero- 
thickness model proposed in Ref. [lij the Neumann BC 



d z 
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(16) 



was used for the potential, corresponding to a vanishing 
electric field at the membrane. Thus the dielectric mis- 
match between the electrolyte and the membrane was 



accounted for only approximatively. When compared to 
the full finite thickness calculation, the agreement was 
poor. To address this issue, a more general Robin-type 
BC was introduced 115] 
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where 



A- m — — d 



0(0+) -0(0"), (17) 



(18) 



is a length scale that contains the membrane thickness d 
and the ratio of the dielectric constant of the electrolyte, 
e, and of the membrane, e m . This BC (with one side 
held at constant potential) was originally developed for 
electrodes sustaining Faradaic current [H, [H, H3, HI] 
or charging capacitively [H, HH]. In that context the 
analog of our membrane is a Stern monolayer of solvent 
molecules or a thin dielectric coating, such as a native 
oxide, on a metallic surface, and A m is denoted As. Note 
that this BC has also been used in Ref. (2lj . 

The modified boundary condition (fT7|) introduces a 
new dimensionless parameter [32| . 



\ m 

K"^ m — — 

Ad 



a 



tF- (19) 



For a blocking or "ideally polarizable" surface, which 
does not pass normal current and only allows capaci- 
tive charging of the double layer, this parameter controls 
the relative importance of the capacitance of the surface 
(here, the membrane) C m = e rn /d compared to that of 
the diffuse part of of the double layer, Cd = e/t. The BC 
(jTTJ) then implies that these capacitances are effectively 
coupled in series in an equivalent-circuit representation 
of the double layer (35|. For a surface sustaining nor- 
mal current, either by electron-transfer reactions at an 
electrode or by ionic flux through a membrane, the sit- 
uation is more complicated. It can be shown that the 
same electrostatic BC (|17p remains valid for a thin di- 
electric layer, as long as it has zero total free charge f34j . 
which is typical for membranes containing a high density 
of fixed counter charge. However, the same parameter 
5 m no longer plays the role of a capacitance ratio. In- 
stead, it controls the effect of diffuse charge on the nor- 
mal current, the so-called "Frumkin correction" to reac- 
tion kinetics in electrochemistry, reviewed in Ref. |36l| . 
Two distinct regimes were first identified in Ref. [321 ] in 
the context of electrolytic cells and recently extended to 
galvanic cells [36|]: (i) the "Helmholtz limit" S m ^> 1, 
where most of the double layer voltage is dropped across 
the surface or membrane and the diffuse-layer has no ef- 
fect on the current, and (ii) the "Gouy-Chapman limit" 
S m <C 1 where the diffuse layer carries all of the voltage 
and thus determines the current. In the Helmholtz limit, 
the Robin BC (fT7|) reduces to the Neumann BC (fl6)l used 
in Ref. [1 51 ] . so that paper analyzed the limit where the 
diffuse charge is small and has little effect on the current. 
In this paper we consider the general case of finite 5 m . 
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III. BASE STATE 

The base state of the problem is a fiat membrane. 
The electric field, assumed to be perfectly aligned in z- 
direction, is then perpendicular to the membrane. In 
the bulk fluid, the system is completely characterized by 
the electrostatic potential 0o( z ) ( or by the field E (z) — 
Eo(z) = —d z (f>o), by the steady-state ion distribution 
Po(z), as well as by the pressure Po(z). Inside the mem- 
brane, an internal electrostatic potential 0™ (z) (and field 
-E™) is present. 

Equation (Tj"2")) is readily solved leading to the charge 
distribution 



P0{Z) = <j KZ 



z > 
z < 



(20) 



Insertion into Eq. (jlip and integrating once yields the 
potential 



<l>o(z) 



I 2 



"fa (* + f) ; *<o - ( } 

Here j m = —j p is the electric current density and 
p(z = 0+) - p(z = 0-) 1 

Pm= g = : 2^J*=° t 22 ) 

represents the jump in the charge density across the 
membrane. We have introduced the notation 



[f]z=a = f( z = a+ ) - fi z = a )> 



(23) 



by which we denote the jump of the field / at position 
z = a. Note that the jump in the charge density p rn at 
the membrane can be interpreted in terms of a surface 
dipole localized on the membrane. The existence of this 
surface dipole is the physical reason for the discontinuity 
of the potential at the membrane. 

At the membrane, the Robin-type BC reads 

A m 9 z 0|2 =o+ = A m <9 z 0| 2=o _ = [4>] z=0 . (24) 

Using Eq. (|21[) . the potential jump at the membrane is 



A 1 2 / j m i 

£>o z =o = — 2 „ 2p„ 



V. 



(25) 



Although the membrane has a zero thickness in this 
model, one can still define an internal field E™ and 
an internal potential <f>™(z). This is done by keeping 
a finite thickness d at first, and then take the limit 
d — > [l5j]. The continuity of the potential at the 
membrane boundaries then implies a constant internal 
field Eq 1 = —[4>o] z =o/d, where the internal potential is 
^(z) = [<t>a]z=oz/d, or explicitly 
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Using Eq. 
density 



24]) . one obtains the jump in the charge 



Pn 



2 v 



-(L + X m ) 



2 + k\ t , 



(27) 



Two remarks on this derivation are in order: first, 
Eq. (|27p illustrates that the asymmetry of the charge dis- 
tribution results either from the accumulation of charges 
due to the applied voltage (a capacitive effect propor- 
tional to V, present also for non-conductive membranes) 
or due to ionic currents across the membrane (a non- 
equilibrium effect proportional to j m present only for con- 
ductive membranes). Second, the expressions for p(z), 
4>(z), p m derived above are independent of the response 
of the ion channels and thus remain unchanged if a non- 
linear ion channel response is used. Only the expression 
for the current j m , that enters as a parameter, will be 
affected by such a nonlincarity. For precisely that rea- 
son, the general form of the base state does not depend 
on the mechanism that has created the current (external 
due to an applied field or internal due to pumps). 

To determine the current density j m at the membrane 
position, we use a linear response approach (for nonlinear 
ionic response, see for instance Refs. (3. l2ll. l39j ) 



e 



-0 ■ 



(28) 



where pP is the chemical potential per particle and G is 
the membrane conductivity per unit surface (across the 
membrane, not in- plane). We assumed equal G for both 
ion species. In the bulk one has j p = — e ® n T d z p p , which 
leads to the usual expression for the (electro-) chemical 
potential, p p = ksT- 1 ^ + e<p. By equating 



G fk B T 



'|z=0 



[po] 2 =o + e[0oUo , (29) 



we finally arrive at the simple expression 
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(30) 



This relation is consistent with the usual electric circuit 
representation of ion channels in a membrane that is sur- 
rounded by an electrolyte 0, EH- As far as the sign 
convention of the currents is concerned, the cathode (to- 
wards where the cations drift) is located at z = —L/2 
and the anode at z — L/2. Thus j m is positive, in accor- 
dance with the usual convention for transport of positive 
charges from the anode to the cathode. Insertion of j m 
into p m yields 



€K 



V 



1 



1 r-, G \ r i 

€K, Z D " 



(1 + j^jyGL) (2 + n\ m ) 



(31) 



The jump in the charge at the membrane is thus posi- 
tive, p m > 0, when A m < £K 2 D/(2G). For reasonable 



5 



parameters for a biological membrane (see section I VI I) , 
p m > holds. A negative charge jump is only possible if 
the electrolyte contains very low salt and in addition the 
membrane is strongly conductive. This is the reason why 
the zero thickness model with the simple BC, Eq. (fT6)) . 
is not realistic for biological membranes, since it corre- 
sponds to the case A m = oo, implying p m < 0. We note 
that by performing the limit A m — > oo in Eq. (|3ip . one 
regains the results for the zero-thickness model of Ref. 
pj| with BC Eq. iflo]). 

To complete the description of the base state, we have 
to consider the total stress tensor 

7-y = -PSij + x) {d lVj + djVi) + e {e.Ej - ~SijE^ (32) 

at the membrane. It contains the pressure, the viscous 
stresses in the fluid and the Maxwell stress due to the 
electrostatic field. We denote by r) the viscosity of the 
electrolyte and by v its velocity field. The electric field 
is given by E = — V(/>. 

In the base state, where the membrane is flat and the 
electric field is oriented in z-direction, from V • r = we 

get d z P Q = %d z ((<9 Z O ) 2 ) = -2p d z <p . By using Eqs. 

(|20l l2Tj) and imposing P(z — > oo) = 0, this is readily 
solved leading to 

P Q ( Z > 0) = -1 ^e-« + ^ e - a "*) , (33) 

and similarly with z — > — z for z < 0. For the stress we 
thus get 

r zz ,o(z > 0) = r zzfi (z < 0) = — ^ j 2 m . (34) 

Note that the stress is constant and is due to the cur- 
rent density j m , with no contributions from the induced 
charges p m . At the membrane, the stress is balanced. In 
addition to the part of the stress tensor due to the elec- 
tric field in the electrolyte, there also is the part due to 
the electric field inside the membrane, already mentioned 
above. For the force balance in the base state however, 
this contribution vanishes is thus not important. 



IV. LEADING ORDER CONTRIBUTION OF 
MEMBRANE FLUCTUATIONS 

In the following, we derive the corrections to the base 
state to first order in the membrane height h(r±). From 
such a calculation, we obtain the growth rate of mem- 
brane fluctuations by imposing the BC for the normal 
stress at the membrane. From this growth rate, we then 
can identify electrostatic and electrokinetic corrections to 
the elastic moduli of the membrane. 



A. Electrostatics to first order in h(r±) 

We use here the quasi-static approach HH by as- 
suming that membrane fluctuations are much slower than 
the characteristic diffusion time tjj = jj^j of the ions to 
diffuse on a Debye length. With the definition of the 
Fourier transform /(kj_,z) = J dr±e^ llc± ' r± f(r±, z), we 
expand the electric field and ion density as 

<f>(k±,z) = <f>o(z) +<l>i(k±,z) , (35) 
p(k±,z) = p (z) + p 1 (k±,z) , (36) 

where kj_ lies in the plane defined by the membrane and 
4>o{z), po(z) are the base state solutions given by Eqs. 
([2U1 |2"T} . The governing equations to this order read 

{dl - k\) <Mk ± ,z) + ^Pi(kx^) = 0, (37) 

{d 2 z -kl) ( Pl (k ± ,z) + | K 2 1 (k ± ,z)) =0, (38) 

and for the particle currents at the membrane one has 
the condition 

DdJ Pl + U 2 ^) = G (—2 [piL=o+ [^}z=o] • (39) 

Eqs. ([351 [Ml) and the BCs at infinity (since L > X D , 
we can use 0i(kj_,±oo) = pi(kj_,±oo) = far from the 
membrane) are satisfied by choosing 

Pi = -T^i (40) 

(and ri^ = —n^ implying pi = en|). It follows that to 
first order in the height, one has a zero flux condition 
at the membrane. Accordingly, the zeroth order solution 
enters in the equations for the first order solution cf>i only 
via the boundary conditions. 

It remains to solve the Poisson equation (|37| . 

(dl - k\ - k 2 ) 0i = . (41) 

We introduce 

I 2 = k\+K 2 , (42) 

with Z _1 as the characteristic length for the electrostatic 
problem near the slightly undulated membrane, and eas- 
ily get (j>i — A T e Tlz for z > and z < respectively. To 
determine the constants A ± , we expand 

9z4>\z=h(T±) = d z 4>Q\ z=h + d z (f>i\ z=h 

= d z <f> \z=o + dz't>o\z=o-H r ±-) 

+d t th\M=o + 0{h 2 ) , (43) 

to prescribe the BC at the membrane in first order. Using 
this expansion in the Robin-type BC, Eq. (|24[) . we get 



0i(kx,2) = 




for z > and a symmetric expression (with e lz ) for z < 0. 
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B. Linear hydrodynamic flow 

When the membrane starts to undulate with small am- 
plitude h(r±), a flow is induced in the surrounding elec- 
trolyte. Using again the quasi-static assumption and low 
Reynolds number, this flow is governed by incompress- 
ibility and the Stokes equation, 



V-v = 0, 
-Vp + rjAv + f = 0, 



(46) 
(47) 



where f is a body force density due to the electric field. 
Introducing the triad [3, H(| of unit vectors (kx,n,t) 
with = k±/k±, n = z and t = ki x n, we get 



d z v z + ikj_ ■ vj_ 
-ik ±P + V (d 2 - k\) v± + f± 
-d z p + V (d 2 -k 2 ± )v z + f z 
V (d 2 - k\) v t + f t 



(48) 
(49) 
(50) 
(51) 



Two forces drive the flow: one is given by the cou- 
pling to the membrane, enters via the BC, and is dis- 
cussed below. The second one is the bulk force due 
to the electric field acting on the charge distribution 
and reads f = QE = —2pV(f). Note that one has 
to use the total charge density, Q — 2p. To lead- 
ing order in the membrane height, this driving force is 
f = — 2poV4>i — 2 / o 1 V0o + 0(h 2 ) with components 

fx = -2po(«)ikx^i(k x ,«), (52) 
f z = -2p (z)d z cl>i(kj_,z)-2pi(k±,z)d z (l>o(z), (53) 



and ft = 0. Because of the latter, Eq. (|5T|) is decoupled 
and trivial. The equations for v z and vi can be decou- 
pled as follows: Using incompressibility, Eq. (|48|) . for the 
perpendicular fluid velocity in Eq. (|49j) . one obtains for 
the pressure 



-rjd z v z 



k f 



rrd z v z 



Insertion into Eq. (|50p yields a single equation for v z 
For z > one has to solve 



(d 2 z -klfv z = 



2 2 / e 2 
-k ± cf)id z [po + -k 0o 

77 V 2 



4 « ra ™e^. (54) 



Imposing v z (z — > ±oo) = 0, the solutions are of the form 



v z = (B^ + C T z) e* k±z + F T e Tlz 



(55) 



for z > and z < respectively. The coefficients 
are determined by the driving force f , but two more BCs 
are needed. At the membrane, continuity of the normal 
velocity imposes 

«,(0+)=«,(0-) = d t h{r ± ) = sh(k x ), (56) 
where we have introduced the growth rate for membrane 
fluctuations, s, from the temporal Fourier representation 
h(t) oc e st . Note that s is also a function of k^. 

The continuity of the tangential velocity, vj_(0 + ) = 
vj_(0 - ) = 0, together with the incompressibility implies 
a second BC for v z , namely 



d z v 



d z v z \ - 



0. 



(57) 



This allows to solve the complete linear hydrodynamics 
problem. With the notation 



eDr)K 4 



(58) 



which quantifies the amplitude of the ICEO flow, the 
velocity and pressure fields in the domain z > read 



v z {z > 0) 
vi(z > 0) 
p(z > 0) 



Mkx) 



s (1 + k±z) e 



h(k±) i kj 



h(k ± ) V 



- hi z 



-k±z 



ak i I k i z 



2k i s — a I k 



fc_L k 2 , z 



a 1 — k I z 



r 



k\z 



I I J I 



ae 



-lz 



lz i ^Pm c -(1 + k)z 

erjl 



(59) 
(60) 
(61) 



The solutions for z < can be obtained by symmetry op- 
erations: v z (z < 0) is obtained by performing the mirror 
operation with respect to the plane defined by the mem- 
brane, z — > — z, in the formula for v z (z > 0). Similarly 



— v±(z < 0) and — p(z < 0) are obtained by doing this 
operation on v±(z > 0) and p(z > 0), respectively. 

In the absence of electric effects, p m = = j m , one gets 
the typical flow induced by a membrane bending mode 
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(4ll . |42| . An additional flow field due to the membrane 
currents arises that has the form of an ICEO flow [l8| . 
A detailed discussion of this effect is postponed to sec- 
tion Efl but we stress that the additional flow is purely 
due to membrane conductivity (since a oc j m oc G), and 
is a non-equilibrium effect. For non-conductive mem- 
branes this flow vanishes but there still is charge accu- 
mulation (p m 7^ 0) in the Debye layers, leading via the 
pressure field, Eq. (|6ip . to corrections to surface tension 
and bending rigidity proportional to p^. 



V. GROWTH RATE OF MEMBRANE 
FLUCTUATIONS 

To discuss the dynamics of the membrane, we still have 
to determine the growth rate s of the membrane. The 
elastic properties of the membrane are described by the 
standard Helfrich free energy 



(62) 



where So is the bare surface tension and Kq the bare 
bending modulus of the membrane. Force balance on 
the membrane implies that the restoring force due to the 
membrane elasticity is equal to the discontinuity of the 
normal-normal component of the stress tensor defined in 
Eq. E2) 



: (0+) + r 2Z (0-) 



dF H 
dh{v±) 



(63) 



We should stress that the coupled electrostatics- 
hydrodynamics problem under investigation can not be 
formulated only in terms of bulk forces, i.e. f and the 
divergence of a stress tensor, because the hydrodynamic 
and Maxwell stress tensors enter the BC (|63|) explicitly. 
For this reason, the force localized on the membrane sur- 
face is a priori unknown, i.e. must be determined by BCs 
for the velocity and the stress. Eq. (|G3[) leads to 

(64) 



[t„,i]m=o = (-£o*i - ^ofcl) fc(ki) : 



which determines the growth rate s = s(kj_) entering the 
stresses. The total normal stress at the membrane is 



T ZZ A 



-P + 2r]d z v z 



\z=h 



(65) 



to linear order in h. Note that we have included here the 
electrostatic contribution stemming from the field inside 
the membrane (with potential 4> m ). This contribution is 
particularly significant in the high salt limit, where effects 
due to the Debye layers become negligible. This inside 
contribution enters with opposite sign than the outside 
contribution because of the difference of orientation of 
the normal (see appendix [Bl and Ref. 19]). 
The outer electrostatic contribution reads 



l(d z M(dM z=0 + hd. 



Po + | (d z <j>of 



where the second contribution vanishes since the term 
in the bracket is a constant (the stress is balanced, see 
above). The electrostatic contribution from inside the 
membrane can be expressed analogously, and the normal- 
normal stress difference at the membrane reads 

[r gz ,i]z=o = -\p]z=o + 2v[dzV z ] z =o 
+e P^oX^OUo - e m [(^<OW)U ±d/2 . (66) 

By means of Eq. ([59]) one easily verifies [d z v z ] z=0 = due 
to the symmetry given above. For the pressure difference, 
Eq. (j6"Tj) implies [p] z =o = 2p(0 + ) and after reexpressing 
I by k± and expanding in powers of k± , one obtains 



[p]z=o = h(ki 



PI, 



Pmjm 

eDn 2 



fc 



€K 



16 



Pmjr. 



(67) 



ere 3 eDn 4 J ■"- L ' " eDn 5 ' 
The electrostatic contribution from the electrolyte reads 

Prnjr 



[(£Wo)(^i)U = fc(k±)8 



eDk 2 



& ) . (OS) 
en 



which exactly cancels the fc_L-independent contribution 
of [p]o- The calculation of the electrostatic contribution 
from inside the membrane is slightly more involved and 
is detailed in Section [5] The result is 



p z c)(^7)L 

-h(k ± )e m (E™y 



=±d/2 
dk 2 , - 



12 



(69) 



where E™ is the zeroth order electric field inside the 
membrane introduced in section [TTT1 

By collecting all the contributions in Eq. (|64[) . the 
growth rate s(kj_) finally has the form 

1 



r]kxs(k x ) = -|(S + AS)fci 



- (K + AK) hi 



(70) 



The electrostatic corrections to the surface tension, 
AE = AE K + AS TO , and to the bending modulus, AK — 
AK K + AK m have been decomposed into outside con- 
tributions (due to the Debye layer, index n) and inside 
contributions (due to the voltage drop at the membrane, 
index m). They are given by 



AE K = _ 4 %-16- Pm> 
en 6 



en 4 D 



AK K 



(71) 
(72) 



for the contributions due to the Debye layers and by 



AE m = - 

AK m = £r, 



m\2. 



z=0 



d 3 



( Tprn \ 2 I _ mi " 
^0 ) 1n Ero 



12 



E%" en" 



(73) 
(74) 
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for the contributions due to the field inside the mem- 
brane. We note that the ratio AK/AH is independent 
of the applied voltage and scales with the square of the 
membrane thickness. 

An independent check of Eqs ([7T1 175)) is provided by a 
direct integration of the lateral pressure profile as shown 
in appendix [Bj This route avoids the consideration of 
hydrodynamics but is limited to the calculation of the 
surface tension correction. 

In Eq. (|70[) . we also obtain a purely non-equilibrium 
correction 



It corresponds to a term proportional to k± in the ef- 
fective free energy of the membrane, which is forbid- 
den for an equilibrium membrane but allowed in non- 
equilibrium. This particular contribution arises due to 
the electroosmotic flows around the membrane as can 
be shown from a simple calculation using the Hclmholtz- 
Smoluchowski equation for the electro-osmotic slip veloc- 
ity near the membrane [l5j |. 

Two particular limits have been considered before: 
first, the case of a non-conductive membrane (J m — 
0, p rn ^ 0) with an arbitrary amount of salt, where we 
recover results obtained recently by Ambjornsson et at 
(l9| . as further discussed in the next section. Second, 
the high salt limit (k — > oo) of a conductive membrane. 
Here, one finds that the correction to the surface ten- 
sion is only due to the inside field, since AS re — > and 
A£ m — > —e m ^f, as has been calculated in Ref. [2j|. The 
corrections to the bending modulus were not considered 
in that reference, and are given from our calculation by 
AK K — > and AK m — > e m -jjd, as in the non-conductive 
case. In this limit, AK m /AT, m = —d 2 /12. 



VI. DISCUSSION 
A. Applications of the model to experiments 

Recently, S. Lecuyer et al. (2f| have investigated a 
pair of nearby membrane bilayers in an electric field by 
neutron reflectivity. The first bilayer was close to the bot- 
tom electrode and used to protect the second one from 
interacting with the wall. Since the bare values of the 
elastic moduli were known from X-ray off-specular ex- 
periments for a similar system [43] (So — 0.5mNm _1 
and Kq ~ 40fcsT), the surface tension correction was 
extracted from the data under the assumption that the 
bending modulus is not affected by the field. The exper- 
iments were performed in an AC electric field at several 
frequencies. For the lowest frequency (10Hz), the elec- 
trostatic correction to the surface tension was obtained 
to be AE ~ -3mNm -1 . 

In the experimentally probed regime of low salt (D2O 
was used as the electrolyte), the electrostatic corrections 



to the elastic moduli depend rather sensitively on both 
the amount of salt and on the ionic conductivity of the 
membrane, as we will see below. Moreover, in the above 
experiment, the correction to the bending modulus was 
not measured. Thus we restrain ourselves to a compari- 
son of orders of magnitude only. The dielectric constants 
are e = 80eo for the electrolyte e m = 2eo for the mem- 
brane. The membrane thickness is typically d = 5nm 
leading to A m — j-d — 200nm. The diffusion coefficient 

of ions is of the order of D = 10 ms, and the viscos- 
ity rj = 10 _3 Pas. The distance between the electrodes 
was about L = 1mm, while the voltage was in the T5V 
range. Assuming that the membrane is non-conductive, 
G = 0, for K = 2 • 10 7 m _1 and V = IV, our model yields 
AE ~ -2-mNm- 1 and AK ~ 190k B T. The model thus 
successfully accounts for the order of magnitude of the 
electrostatic correction to the surface tension. However, 
it also shows that the bending modulus increases about 
five times. In order to obtain an experimental test of 
the model, it would be interesting to measure the correc- 
tion to the surface tension and to the bending modulus 
simultaneously. We would also like suggest to carry out 
experiments in which the applied electric field or the ionic 
strength would be varied. Another interesting possibility 
would be to study membranes of different conductivities 
or thicknesses in an applied electric field. 

A second field of application of the model are active 
membranes, which are (artificial) lipid vesicles containing 
ionic pumps such as bacteriorhodopsin @, H, • In these 
experiments, no external electric field is applied. Instead 
the pumps are activated by light to transport protons 
across the membrane. In Ref. [9|, a lowering of the mem- 
brane tension produced by the activity of the pumps has 
been reported, which could be due to an accumulation of 
charges near the membrane, as discussed here. The speci- 
ficity of that experiment is that this charge accumulation 
would result from the activity of the pumps rather than 
from an applied electric field. However, it is difficult to 
make a precise comparison between the experiments and 
the present theory, because only the correction to the 
surface tension is accurately measured and many aspects 
of the transport of ions are unknown. Nevertheless, if we 
assume that the passive state of that experiment corre- 
sponds to a non-conductive membrane (G = 0) and the 
active state to a membrane with G = 105m -2 , and if 
we use a typical transmembrane potential of the order of 
50mV, we get the same order of magnitude for the ob- 
served tension lowering, 3 • 1Q — 'Nm -1 , if we account for 
the rather high amount of salt with k ~ 5 ■ 10 8 m _1 . We 
also find that there is no measurable difference for the 
bending modulus between the active and passive state, 
as observed experimentally. The model predicts that a 
current density of j m ~ lAm~ 2 arises when the pumps 
are active, which corresponds to an overall current of 
lpA on a vesicle of size l/zm. To better compare to the 
model, again it would be desirable to have experiments 
in varying conditions (ionic strength and conductivity of 
the membrane, for instance). Another interesting pos- 
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FIG. 1: Panel a) shows the electrostatic corrections to the 
surface tension (in units of N-m -1 ) and panel b) those to the 
bending modulus (units J) as a function of k (units m _1 ) 
in the non-conductive case, G = 0. Dashed lines: contribu- 
tions due to the Debye layer (AE K and AK K respectively). 
Dash-dotted lines: contributions due to the field inside the 
membrane (AE m and AK m respectively). Solid lines: sum of 
both corrections. The figure was made with the parameters 
given in the text and V = IV, L = 1mm. 

sibility would be to measure the membrane current and 
the transmembrane potential in the course of the exper- 
iment, for instance using patch-clamp techniques. 

B. Effect of salt and membrane conductivity 

With the expressions for the jump of the charge density 
at the membrane, Eq. (|3Tj) . and for the current density, 
Eq. (|30p , we can discuss the corrections to the membrane 
elastic constants given by Eqs. (|7T|) - (|71)) . In particular 
we obtain the dependance of these elastic moduli on the 
ionic strength of the electrolyte, and on the ion conduc- 
tivity of the membrane G. Fig. []J,) displays separately 
the contributions to the surface tension as a function of 
the inverse Debye length k. The dashed line represents 
the contribution from the Debye layers, the dash-dotted 
line represents the contribution from the field inside the 
membrane and the solid line is the sum of both con- 
tributions. The value of the inverse Debye length re, 
varies from re ~ 10 6 m _1 (Xjj = lpm) for pure water 
to re — 3.3 • lOW 1 (Ajj ~ 0.3nm) for 1M NaCl @. Fig. 
HJd) shows the respective contributions to the membrane 
bending modulus. In this figure we have assumed that 
the membrane is non-conductive (G = 0). As shown in 
Fig. [TJ the contributions from the Debye layers dominate 
for low salt (re < 5 • 10 6 m _1 for AS and re < 10 8 m _1 for 
AK). For high salt, the contributions from the mem- 
brane dominate, and approach the limiting values dis- 
cussed above. In the case of zero conductivity, both the 
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FIG. 2: Electrostatic corrections to the surface tension (panel 
a) and to the bending modulus (panel b) as a function of ft 
in a slightly conductive case (G = O.lSm" 2 ). Dashed lines: 
contributions due to the Debye layer. Dash-dotted lines: con- 
tributions due to the field inside the membrane. Solid line: 
sum of both corrections. Parameters as in previous figure 
except for G. 



Debye and the inside contribution to the surface tension 
are always negative, and there is good agreement with 
calculations for non-conductive membranes (l5l . [l^ | . Fig. 
[2] displays the corrections to the elastic coefficients in 
the case where a finite membrane current j m is present, 
induced by a small conductivity G = O.lSm -2 , for oth- 
erwise unchanged parameters. We find that this rather 
small conductivity has already a large effect on both 
moduli: first, the effect of the Debye layers on the surface 
tension is suppressed and the contribution from the mem- 
brane is dominating. Second, the overall contribution 
gets relevant for higher salt than in the non-conductive 
case. The effect on the bending modulus is even more 
significant: although the Debye contribution is still dom- 
inating for about re < 10 m , it is much smaller in am- 
plitude than in the non-conductive case (around 10 -20 J 
compared to 10~ 17 J at n ~ 5 • 10 6 m _1 ) and furthermore 
becomes non- monotonous fl5j . 

The effect of membrane conductivity is highlighted in 
Fig. [3j where the total contributions to surface tension 
(panel a) and bending modulus (panel b) are shown as 
a function of the inverse Debye length re for conductiv- 
ities in the range G = 0.01 — lOSm -2 , and otherwise 
unchanged parameters. We find that, except for the high 
salt limit, the bending modulus correction tends to be 
reduced by increasing membrane conductivity. To give 
some numbers, for a non-conductive membrane (G = 0) 
and V — IV, the jump in the charge density for re = 
2 • 10 7 m _1 is p m = 1.5 ■ 10 23 ^. Already a small value 
of the conductivity as G — 0.15m -2 halves the charge 
density to p m = 8.6 • 10 22 ^- and creates the current den- 
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FIG. 3: Electrostatic corrections to the surface tension (panel 
a) and to the bending modulus (panel b) as a function of « for 
different membrane conductivities. Parameters are as in the 
two previous figures except for G: solid line G — 0; dashed 
line G = 0.01; dash-dotted line G = 0.1; dash-two-dots line 
G = 1; dotted line G = 10 in units of Sm~ 2 . 



sity j m = 3.7 • 10 17 -Sr, or 60 • 10~ 3 A- Conductivities 
of biological membranes can be as high as G = lOSm -2 , 
which is the value for a squid axon, corresponding to a 
density of potassium channels of 0.5/xm -2 [4(. 

We also mention that the distance between the elec- 
trodes L, i.e. the confinement, is a relevant parameter 
and influences the shape of Figs. [T]and[2l Here we have 
used a macroscopic distance (L = 1mm), corresponding 
to the experiments mentioned above. If L was instead 
of the order of microns the suppression of the bending 
modulus correction due to conductivity would be much 
less pronounced and the corresponding figure would be- 
come similar to the one given in Ref. [l5(. Moreover, 
for high enough membrane conductivity, the Debye layer 
contribution to the surface tension can become positive, 
i.e. stabilizing. In fact, the sign of AS K is governed by 
a factor — r^r ■> ^^T-> ■ In a way similar as discussed 

(eDK 2 +2GL) J 

in section ITTT1 concerning the sign of p m , for G > £ 2 ^ K 
there is a sign change, rendering the correction positive 
for small k. However, the denominator containing the 
distance L between the electrodes suppresses this effect 
for macroscopic distances. It can be seen only if L is 
small, e.g. L = 1/zm as used in Ref. [IJ[. We note that 
the micron scale is particularly relevant to experiments 
with cell membranes submitted to electric fields [22|. It 
is also relevant to experiments that one could propose to 
test these ideas using microfluidics devices. 



C. Membrane instability 

Since the corrections to the membrane surface ten- 
sion are typically negative (with the exception mentioned 
above), they can overcome the bare surface tension So. 
At this point, an instability towards membrane undula- 
tions sets in [23|]. Our theory is able to go beyond pre- 
vious modeling of this instability (still for early stages of 
the instability) , which were limited to the high salt limit 
and did not include electrostatic corrections to the bend- 
ing modulus or hydrodynamic effects associated with the 
modulus r. The linear growth rate of the membrane 
fluctuations is given by Eq. (|70|) and has the form 

Vs(k) = -^ eff k-T K k 2 -^K eff k 3 , (76) 

where we have written simply k for k± and introduced the 
effective surface tension and modulus, S e /j = So + AS, 
K e ff = K + AK. Fig. [3^) shows this growth rate, or 
dispersion relation, in rescaled units where we scaled the 
wave vector by k, k' — k/n, and the time by the typi- 
cal time for ions to diffuse a Debye length, td = 7573- 
The parameters are the same as in the previous sections 
for a non-conductive membrane, i.e. G = 0. The con- 
trol parameter is the external voltage V. Fig. [4ji) shows 
the growth rate for three different levels of the voltage: 
the dashed line is for V = 0.7V, which lies below the 
threshold of the instability, all wave numbers are damped 
and the membrane is stable. The solid and the dash- 
dotted line correspond to V = 0.75V and V = 0.8V and 
are above threshold. A certain window of wave numbers 
k g]0, k max (V)[ has positive growth rates and the mem- 
brane is thus unstable. This window gets larger with 
increasing voltage. The linear growth will be dominated 
by the maximum of the growth rate defining the fastest 
growing wave number kf g . Given Eq. (|76|). one easily 
calculates 

*/. = H (-r« + ^|%%) , (77) 

for S e // < 0. 

The same information given by the dispersion relation 
can be expressed by the so-called neutral curve which 
is shown in Fig. [3)d). This curve, given by the solid 
line, separates the negative (below) from the positive 
(above) growth rates in the control parameter- wave num- 
ber plane. If the voltage is below the section of the 
neutral curve with the voltage-axis, the system is stable. 
Otherwise a certain band of wave numbers is unstable. 
The position of the fastest growing mode kf g is given by 
the dashed line. 

Since we have the dispersion relation in analytical 
form, in principle one has formulas for all relevant observ- 
ables like the threshold voltage V c . In terms of the system 
parameters, however, they are quite lengthy. The thresh- 
old voltage is given by the change of sign of the leading 
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FIG. 4: a) The renormalized growth rate or dispersion rela- 
tion, tds, as a function of the rescaled wave number k' = 
k±/n for three voltages: V = 0.7V (dashed line), V = 0.75V 
(solid line), V = 0.8V (dash-dotted line), b) The neutral 
curve (solid line) separating the regions of s > and s < 0, 
and the fastest growing wave number kj g (dashed line) in 
the plane voltage vs. rescaled wave number k' — k±/n. 
Parameters as previously except: no conductivity, G = 0; 
k = 2 ■ 10 7 m _1 ; E = lmNm -1 ; K = Wk B T. 



order contribution in s(k±) 
it has the simple form 



In the non-conductive case 



V 2 (G = 0) = 



S rf(2 + n\ m ) 2 



(78) 



Since both p m (and j m in case of G ^ 0) are proportional 
to the voltage, as expected the critical voltage scales like 
V c oc ySo- I n the limit of small membrane conductivity, 
one gets to leading order (using that L is macroscopic) 



V 2 = V 2 (G = 0)[l + 



4GL 
eDn 2 



(79) 



Thus membrane conductivity increases the voltage value 
needed to cross the instability. In the limit of high salt, 
K — > oo, one regains the known result V 2 — Ylod/e m . The 
typical wavelength of the membrane undulations above 
threshold (i.e. the one of the fastest growing mode) for 
parameters as in Fig. U]is of order A = ~ 0.25/^m ~ 
12.5Ad, so several times the Debye length. 



D. ICEO flows 

We now discuss the form of the fluid flows which arise 
near the membrane when it is driven by ionic currents. 
Fig. shows the flow field for a high membrane conduc- 
tivity and low salt, in the regime where the membrane is 
unstable due to the electrostatic correction to the surface 
tension and thus starts to undulate. This figure was gen- 
erated by selecting the fastest growing wave number fc/ ff , 
defined in section IVIC1 and using the respective max- 
imum growth rate s(kf g ). Since this wave number has 
the fastest growth rate in the linear regime, it will dom- 
inate the initial behavior. The shape of the membrane 
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FIG. 5: Representation of the flows around the membrane 
beyond the instability threshold. The orientation of the elec- 
tric field is towards negative values of z. Panel a) shows the 
flow generated by the membrane bending instability (terms 
proportional to s in Eqs. (J59J [60j ) - Panel b) shows the ICEO 
flow (terms proportional to a in Eqs. (|59l [60}). Finally, 
panel c) shows the actual flow, which is the superposition of 
the former two and results in a strong flow near the mem- 
brane, oriented parallel to the surface. Both axes are scaled 
by the Debye length Parameters are as in previous fig- 

ures except V = 3.165V, « = 10 7 m _1 , G = lOSm -2 and 
L — 10pm. 



undulation is represented as the black solid curves in all 
plots of Fig. [5j The resulting flow, shown in Fig. Et) 
is a superposition of two distinct flows: first, the typical 
flow associated to a membrane bending mode [U 52] as 
shown in Fig. [5ji). This contribution corresponds to the 
terms proportional to the growth rate s in Eqs. (|591 I60[) . 
Second, the flow associated with the remaining terms in 
Eqs. ([591 lrj(J|) , proportional to a. This contribution yields 
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the typical counter-rotating vortices of an ICEO flow [l8| , 
as shown in Fig. [SJa) . Clearly, the superposition of these 
two flow contributions, as shown in [5b), results in a par- 
allel flow close to the membrane, in contrast to the usual 
bending mode flow given by Fig. 

Since the jump of the charge density p m > for biolog- 
ical membranes (and j m > by definition), the induced 
flow occurs for this case in the same direction as in stan- 
dard ICEO flows. Note that an inverse ICEO flow was 
obtained in [lj|, due to the opposite sign of p m obtained 
with the simple but unrealistic BC Eq. (fTB")) . Also, the 
situation of ICEO flows is less general as suggested ear- 
lier: For most parameters (modest conductivities, not too 
low salt) the flow generated by membrane bending is usu- 
ally dominating and hides the small ICEO contribution. 
This is due to the fact that the former is proportional 
to s which has contributions cx p^, while the ICEO flow 
is oc pmjra <C p^. Thus, to see the situation given by 
Fig. [5l a high membrane conductivity G is needed. Sec- 
ond, one needs low salt, since otherwise the membrane 
instability is shifted to very high voltages. Also, since 
for macroscopic electrode distances L (of order millime- 
ter) and high conductivity the voltage needed to induce 
the instability is very high, we have used a microscopic 
electrode distance L — 10pm. While it might still be pos- 
sible to see these flows for higher salt and macroscopic 
electrode separations, such situations will be clearly far 
beyond the Debye-Hiickel approximation used so far. 

The ICEO flows near the membrane could also become 
relevant once the system has reached a steady-state. In- 
deed in the case of lipid vesicles for instance, non-linear 
effects associated with the conservation of the number 
of lipids on the vesicle [23| guarantee a saturation of 
the membrane fluctuations (for not too high voltages 
that might lead to vesicle rupture), as compared to the 
case of the planar membrane considered here. Since the 
membrane fluctuations are confined by non-linear effects 
and become quasi-stationary in the long time limit, the 
system can reach a well defined non-equilibrium steady 
state. In this non-equilibrium steady state it might well 
be possible that fluid flows still arise due to ionic cur- 
rents going through the membrane, while the initial flow 
associated with the membrane bending mode has disap- 
peared. 



VII. CONCLUSIONS AND PERSPECTIVE 

This paper offers a route to describe capacitive effects 
near a conductive lipid membrane while keeping the sim- 
plicity of the zero thickness approximation on which most 
of the literature on lipid membranes is based. These ca- 
pacitive effects are the main player in the corrections to 
the elastic moduli of membranes driven by an electric 
field or by internal pumps or ionic channels. The present 
theory goes beyond available descriptions by including 
non-equilibrium effects which arise due to ionic mem- 
brane currents. These ionic currents have a similar form 



as the ICEO flows studied in the context of microflu- 
idics and can modify the fluid flows around the mem- 
brane from usual bending dominated flow towards flow 
concentration close to and parallel to the membrane. 

Our approach is sufficiently simple to be the starting 
point for further generalizations, which could include var- 
ious nonlinear effects: non-linear elastic terms associated 
with the membrane or the cytoskeleton in case of a bio- 
logical membrane, non-linear current- voltage response of 
the channels. Also density fluctuations of the ion chan- 
nels and the experimentally simpler case of an AC electric 
field should be investigated. Further theoretical work is 
also needed to extend the model to higher voltage where 
the Debye-Hiickel approximation breaks down. In fact, 
this approximation only holds when the potential satis- 
fies everywhere the condition J.^ T <C 1 [16]. For the 
experiments on supported membranes discussed above, 
one finds that the charge accumulation on the membrane 
is too large for this approximation to hold, since one has 
jzrm — 2. It would thus be relevant to solve the nonlin- 
ear Poisson-Nernst-Planck equations and otherwise pro- 
ceed similarly as proposed in this work, in order to de- 
scribe the behavior of membranes surrounded by high 
charge densities. Furthermore, it would be interesting to 
investigate a model suitable for small system sizes, since 
much of the results of this paper are based on the as- 
sumption that the system size L is much larger than all 
other length scales in the problem, as well as for more 
realistic boundary conditions at the electrodes. 
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APPENDIX A: CONTRIBUTION OF THE 
INTERNAL FIELD TO THE MAXWELL 
STRESSES 

Here we are interested in the contribution of the inter- 
nal field inside the membrane to the stress, which is 

[t£i]*=o = e m P z C)(^r)L= ±d /2 ■ (Al) 

Since the internal field at zeroth order is constant inside 
the membrane due to the symmetry of the problem, this 
expression simplifies into 

[rZ,ih=o = -2e m E™ (d z ^) ]z=+d/2 . (A2) 

The first-order field in the membrane is given by (use the 
symmetry or cf. Ref. [l5[ for details) 

<l>?(k ± ,z) = W(k ± ,d/2)-£^ (e fc ^ + e~ k ^) ,(A3) 
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which leads us to 



r zz,l\z=0 



-2e m E^[k±,-)k 



£ k ± d _ ^ 



(A4) 



To obtain an expression for <f>™ (fcx, f), to linear order 
in h, one can write [Til ] 

^(fci. = <f>i(k±, - ft - W| z=(i/2 ■ (A5) 



For the outside potential at the membrane we can ap- 
proximately use <f>i(k±,z) — —~p m here (since the 
exponential decay like e~' z starts at the membrane). Af- 
ter expansion in k± (and assuming L > d) we get 



h{k±) 



rpm i P m l 2 



(A6) 



Using Eq. (|A4|) . for the stress difference, [r™ 1 ] z =o, up 
to order /c^ this exactly yields Eq. (|69|) given above. 



APPENDIX B: ELECTROSTATIC 
CONTRIBUTION TO THE SURFACE TENSION 
FROM THE STRESS TENSOR 

In this appendix, we give an alternative approach for 
the derivation of the membrane tension, which avoids 
solving for the fluid flow around the membrane. Here 
the tension is expressed as an integral over the lateral 
pressure profile deviation or, equivalently, over the excess 
lateral stress [44| , 

Let us call S a closed surface englobing the membrane 
with the normal vector field n. We choose x to represent 
the direction of the lateral stress. The force acting on 
the surface S in the x-direction can be calculated from 
the stress tensor defined in Eq. ([32)) as 



Fr 



x • r ■ n dS. 



(Bl) 



Since r is divergence free, the surface S can be deformed 
into an arbitrary other surface englobing 5, for conve- 
nience to a cube of size L. It is easy to see that the 
integral in Eq. (|B1|) is non-zero only on the faces of the 
cube with the normal along ±x. With dS = Ldz and 
AE = F x /L on the face where n = +x, we arrive at fl9j 



r-L/2 

AE = / r xx (z)dz. 

-L/2 



(B2) 



Eq. (|3"2"j) implies t xx (z) = -P (z) — ^(d z 4> ) 2 , where 
4>o( z ) is the potential and Po(z) the pressure profile in 
the base state. The resulting tension AS is identical 
to Eqs. (41-43) of Ref. fisj ] . where it was expressed as 
a sum of two terms, Eo and Ei. A cancellation of the 
dependance on L occurred in the sum of these two terms, 
as it should be. The present derivation fully justifies 
this point since the choice of the deformed surface was 
arbitrary by construction. Using Eqs. (|2T] [33j IB2j) . we 
obtain 



AE K = -A^% - 16 



Prnjrn 



(B3) 



which is exactly the result for the contribution of the 
Debye layers to the tension using the hydrodynamic ap- 
proach, Eq. (|7ip . A similar calculation gives contribution 
to the tension from inside the membrane, which reads 



AE r 



(B4) 



Note that the same method could be used to calculate the 
change of spontaneous curvature induced by an electric 
potential in the asymmetric case [l9j . 
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